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^ ! Abstract 

[ We study the fluctuations of the stress tensor for a massless scalar field in two and 

> ■ 

\Q ' four-dimensional Minkowski spacetime in the vacuum state. Covariant expressions for the 

(N : 

stress tensor correlation function are obtained as sums of derivatives of a scalar function. 
These expressions allow one to express spacetime averages of the correlation function as 
finite integrals. We also study the correlation between measurements of the energy density 

O 

along a worldline. We find that these measurements may be either positively correlated 
or anticorrelated. The anticorrelated measurements can be interpreted as telling us that if 

£> ! 

. !_h ' one measurement yields one sign for the averaged energy density, a successive measurement 

x : 

with a suitable time delay is likely to yield a result with the opposite sign. 
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I. INTRODUCTION 



Because physically realizable states in quantum field theory are not eigenstates of 
the stress tensor operator, quantum stress tensor fluctuations are a universal feature of 
quantum fields. These fluctuations can have physical effects, including Casimir force 
fluctuations 0, 0, H, EL radiation pressure fluctuations E|, and pas sive fluctuations of 
the gravitational field 0, 0, i, ©, Efl EH, El El, El El, EMl, El M HI IH HI HI . 
Passive fluctuations of gravity are those driven by fluctuations of the matter field 
stress tensor, as opposed to the active fluctuations due to the quantum nature of 
gravity itself. The quantum stress tensor correlation function is singular in the limit 
of coincident points. However, this does not prevent us from obtaining physically 
meaningful results for observable quantities, such as the luminosity fluctuations of 
a distant source seen through the fluctuating spacetime [2£j. These observables are 
expressed as spacetime integrals of the correlation function, which can be defined by 
an integration by parts procedure. Alternatively, one could use other approaches, 
such as dimensional regularization p5| . 

In general, the stress tensor correlation function can be decomposed into three 
terms: a "fully normal ordered" term which is state dependent, but free of singular- 
ities, a vacuum term which is singular, but state independent, and a "cross term" 
which is both singular and state dependent. In many situations, one is interested in 
state dependent effects, so the vacuum term can be ignored. For example, radiation 
pressure fluctuations in a coherent state arise solely from the cross term . However, 
this does not mean that the vacuum term is devoid of any physical content. 

The main purpose of this paper is the search for such content. Here we will be 
concerned with a free, massless scalar field in Minkowski spacetime, and its stress 
tensor correlation function in the Minkowski vacuum state. In a previous paper pfil ]. 
we studied the subtle stress tensor correlations in non- vacuum states created by mov- 
ing mirrors in two-dimensional flat spacetime. One of the key results of the present 
paper will be the derivation of a covariant expression for the correlatio n fun ction as 
a sum of total derivative terms. This expression will be given in Sect. Ill Al for two 
dimensions and in Sect. IIIlAl for four dimensions, with the details of the derivations 
presented in Appendices |X] and |Bj respectively. We wi ll d iscuss spacetime averages 
of the energy density correlation function in Sects. iHBl and UlI Bl and averages along 
a worldline in Sects. Ill CI and IIII CI The results will be summarized and discussed in 
Sect. IIVI Units in which h = c = 1, and a spacelike metric signature will be used 
throughout this paper. 



II. TWO DIMENSIONS 

A. Covariant Stress Tensor Correlation Function 

We will be concerned with the stress tensor correlation function 

C^apix,^) = (: T^(x) : : T aP (x') :) (1) 
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for a massless, minimally coupled scalar field in two-dimensional Minkowski spacetime 
in the vacuum state. Here : T^ u (x) : is the normal ordered stress tensor operator, so 
(: Tfj, u (x) :) = 0. We especially seek an expression for C lxva p{x, x') as a sum of terms, 
each of which is a total derivative of a function with at most logarithmic singularities 
as x' — > x. This will allow us to define integrals of the correlation function by 
integration by parts. 

Such a form is derived in Appendix |XJ where it is shown that 

C^ af3 (x, x') = — -8 df, d u d a dp fx -2 g aj3 DDf 2 

do4 7T z L 

+ {Qua 9u/3 + fiW? gva) 00 f 2 + 2 (g a p d v + g^ d a dp)Uf 2 

- {g av dp. dp + g afl d v dp + g Pv d a + g^ d v d a ) Df 2 , (2) 

where 

h = ln(Ax 2 /£ 2 ) , (3) 

and 

h = hi 2 (A* 2 /£ 2 ) , (4) 

where t is an arbitrary constant with dimensions of length. The correlation function 
is independent of the choice of £. Here □ = <9 M <9 M is the wave operator, and Ax 2 = 
(x M — — x'fj). Because <9 M fi = dfi/dx 11 = —d^ f\ = dfi/dx" 1 , the correlation 

function, Eq. (J2J), can be written in several equivalent forms. 
The energy density correlation function becomes 

C(x, x') = C tm , = Of fx = d\dl h ■ (5) 

Note that none of the f 2 terms contribute in this case. This expression allows us 
to compute the mean squared average energy density. Let g(t) be a time sampling 
function, and h(x) be a spatial sampling function. Then we define the averaged 
energy density operator as 

p = J dtg{t) J dxh{x) :T tt : . (6) 

The mean square of this operator is 

C= (p 2 ) = [ dtg{t) [ dxh(x) [ dt'g{t') [ dx' h{x') C{x,x') . (7) 



If we insert Eq. (JSJ) into the above expression, and then integrate by parts, we can 
write 

C* = ---L j dtg{t) Jdt'g\t') Jdxh(x) J dx' h(x') h . (8) 

In the limit that the width of the spatial sampling function goes to zero, h(x) — > S(x) 
and we obtain 

C = - JL / dfg(t) I dt'g(t') \n[(At) 2 /£ 2 ] . (9) 
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B. Averaging over Space and Time - 2D 



Rather than using Eq. Q, in some cases we can also directly evaluate the integral 
in Eq. (JJj) using contour integration methods. For the explicit examples to be treated 
in this paper, the latter approach is more convenient. The energy density correlation 
function, Eq. (jSJ), can be expressed as 

C(x x') - (^ 2 + Ax 2 ) 2 + 4AtW 

CM ~ 4^ (At 2 - Ax 2 ) 4 ' (10) 

where At = t — t', and Ax = x — x' . In this subsection we will sample this correlation 
function in both space and time with Lorentzian functions of width a in t and t', 
and (3 in x and x'. Further, let the spatial sampling functions coincide, but let the 
temporal ones be displaced by t . 
Let 

/oo poo poo poo 

dtg L (a,t+t ) / dt' g L (a,t') / dxg L ([3,x) / dx' g L (/3, x 1 ) C(x, x') 
-oo J —oo J —oo J —oo 

(11) 



where 



and 



-oo 

Now let t —>■ t — t , so that we have 



g L {a,t)= a , (12) 
ir{t z + or) 

dtg L (a,t) = l. (13) 



C(t ) = / dtg L (a,t) / dt' g L (a,t') / dxg L ((3,x) 

J —oo J —oo J —oo 

/oo 
dx'g L (P,x') C(t — t' — t ,x — x') 
■oo 

/oo poo 
dr g L (a,r) / dp g L (b, p) C{r - t , p) , (14) 
-oo J —oo 

where a = 2a, b = 2(3, r = t — t' and p = x — x'. In the last step, we used the identity 

/oo /*oo /*oo 

dtg L (a,t) j dt' g L {a,t')F(t-t') = J dr g L (a,r) F(t) . (15) 
-oo J —oo J —oo 

We may do the integral on p first, by contour integration. The integrand has 
simple poles at p = ±ib and fourth order poles at p = ±(r — t ). We choose a contour 
in the upper half-plane which avoids the fourth order poles, the contour C\ in Fig. ^ 
In fact, we could use other contours such as C2 and still obtain the same answer. Even 
if we chose a contour which enclosed either of the fourth order poles, our answer for 
the real part of the integral would still be the same. This is because the contribution 
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of either of these poles, the result of integrating around the closed circular paths, is 
pure imaginary. Note that the straight segments and the semicircular segments of 
C\ each contain real terms which diverge as the radii of the semicircles go to zero. 
However, these terms cancel when the straight and semicircular contributions are 
added. The divergent terms on the straight segments are the boundary terms that 
would arise from integrating by parts along these segments only. Thus integration by 
parts along the straight segments and discarding the boundary terms produces the 
same result as integration along the complete contour. 
In any case, using the residue theorem we obtain 

A ( h \rt f \ [(r - t ) 2 - b*f - 4b* (r-t ) 2 

dp 9L{b,p) C(t -t ,p) = 3 . 16 

4vr 2 [(r - t ) 2 + 6 2 ] 

The subsequent r-integration was performed and yields 

A/, s [t (t + 2a + 2b)-(a + b) 2 ][t (t -2a-2b)-(a + b) 2 ] 

^\to) — ri • U<J 

4vr 2 [t 2 + (a + &) 2 ] 

(This and several other calculations in this paper were done using the public domain 
algebraic manipulation program MAXIMA.) In the special case when t = 0, we 
simply have 

d{0) = ^TW (18) 

Let us define 

K(t ,a,b) = ^-. (19) 
In general, we have that C(to) = C(—t ). From Eqs. (fTTj) and (fl9*)l. we find that 

/oo 
K(to,a,b)dto = 0, (20) 
-oo 



and similarly 



K(t ,a,b)dt = 0. (21) 



This result tells us that positively correlated regions (K > 0), and anticorrelated 
regions (K < 0) have equal weight. 



C. Sampling along a Worldline - 2D 

In this subsection, we shall specialize to the case of sampling along a worldline, 
i.e., we will effectively set the width of the spatial sampling function to zero. Define 
a normal-ordered smeared stress tensor operator by 

/oo 
dtg(t,t ) :T tt (t) :, (22) 
-oo 
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FIG. 1: Some possible integration contours for Eq. Q16JI are illustrated. There are two 
simple poles on the imaginary axis, and two higher order poles on the real axis. Both types 
of poles are denoted by the letter X. The contours C\ and C% both yield the same result 
for the integral. Integration around either of the poles on the real axis (dashed line circles) 
give an imaginary result, so the real part of the integral is independent of whether these 
poles are enclosed or not. 

where g(t,to) is a sampling function whose peak is at t = to- Although (S) = in 
the vacuum state, (S 2 ) ^ 0. From Eq. 0, we have that 



(S 2 ) = dt dt'g(t,t )g(t',t )C(t,t') 

J —oo J —oo 

I poo POO 

= ~W1 dtg(t,t ) dt'g(t',to)\n[At 2 /£ 2 }. (23) 

•J —OO %J —oo 

The case we want to consider is two regions of time-sampled energy density which are 
allowed to initially coincide but which are then gradually separated from one another. 
One sampling function has its peak at t' = and the other at t — to. We want to 
imagine sliding these regions away from one another (see Fig. |2J), and examine the 
behavior of the vacuum correlation function as we vary to- 
With Eq. (J22), we can write 

roc pec 

(S(t o )S(0)) = dt dt'g(t,t )g(t',0)C(t,t') (24) 
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increasing to 




FIG. 2: Two sampling functions which initially coincide and are then gradually separated. 

/oo /*oo 
dt dt' g(t,t )g(t',0)\n[At 2 /£ 2 ]. (25) 
-oo J — oo 
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This represents the smeared energy density correlation function for two displaced 
regions along a worldline. We can normalize this quantity by defining 

(S(t o )S(0)) 

K(t ) = {S2{0)) . (26) 

As an example, we take the sampling function to be a Lorentzian. If we set b = 
and a = 1 in Eq. (fT7jl . we find 



K(t , (tg - 2t - l)(tg + 2f - 1) (1 - K + ig; 

= W^f = + ' (2?) 



The choice of b = corresponds to sampling in time only, with displaced sampling 
functions. A plot of this function appears in Fig.|Sfa). The plot is somewhat deceiving 
because there is actually a second positive peak which, on the scale of the plot, is 
too small to be seen. However, it must be there since C{t ) ~ 1/(4 vr 2 ^), as t — ► oo, 
and hence K(to) has to approach from above for large to- The magnified view in 
Fig. [3fb) reveals the second positive peak. We can also see this by computing the 
extrema of K(to) using 



4to( tg-10tg + 5) 
(tl + if 



K%) = _ ^-m-toj (2g) 
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K(t) 




FIG. 3: The graph of K (to) versus to for a Lorentzian sampling function, in units with 
a = 1. Here we have chosen b = 0, so the sampling is in time only. Here (a) shows the 
overall form of K(to), but on a scale which does not reveal the final maximum. This peak 
is revealed on a smaller scale graph, (b). 

One finds that K'(t ) = at: t = (first maximum), t pa 0.73 (minimum), and 
to pa 3.1 (second maximum). 

As a second example, consider a compactly supported sampling function of width 
a with g = g = 0&tt = to±a/2. A simple choice of function which has this form is 

30 

g(t, t ) = g(t - to) = -g (t - to - a/2) 2 (t-t + a/2) 2 . (29) 

a 

The second derivative of this function is 

30[12 (i io f-^ 

and 

(S 2 (0)> = ^_. (31) 
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K(to) 
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FIG. 4: The graph of K{to) versus io f° r the compactly supported sampling function given 
by Eq. (|29|). in units where a = 1. 

Using Eqs. (J2HJ) . (J2EJ), (pH2j), and (j3~TJl . one can evaluate -K'(to), which is plotted as 
a function of to m Fig- HI Note that the number of maxima and minima of K(to) 
for the compactly supported sampling function, given in Eq. (J2HJ), is the same as 
for the Lorentzian sampling function shown earlier. However, for the compactly 
supported sampling function case, the second maximum is much more pronounced. 
A calculation also shows that for both the Lorentzian and the compactly supported 
sampling functions, we have that 



We will show in Appendix that this is true for arbitrary smooth sampling func- 
tions. In this appendix, we also prove that 



This establishes that the behavior illustrated in Fig. H] is independent of the details 
of the sampling function. The fact that (>S 2 (0)) > implies that nearly overlapping 
regions are positively correlated with one another. As to increases, the correlation 
is replaced by anticorrelation, as shown by the negative minimum in K(t ). This 
anticorrelation implies that if we measure positive energy in a given region, there 
must be negative energy found nearby. Finally, when the regions are sufficiently 
separated, the positive correlation returns, as evidenced by the final positive peak in 




(32) 



(S 2 (0)) > 0. 



(33) 
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Fig. |3J One can understand why disjoint regions must be positively correlated from 
the fact that C(x,x') > 0. When x ^ x' everywhere in the range of integration, then 
the integral for C is well defined as an ordinary integral, and must be positive. On 
the other hand, when we must integrate through points where x = x', then C(x,x') 
becomes defined only as a distribution, and the integration by parts procedure can 
produce a negative result. 



III. FOUR DIMENSIONS 



A. Covariant Stress Tensor Correlation Function 

In this section, we consider the vacuum stress tensor correlation function in four 
dimensions. The covariant form of this function is derived in Appendix [Bl with the 
result 

C^ vaP {x,x') = - r ..] n j \^d^d u d a dpUUf 2 + Qg^g al3 U i f 2 

0144U7T* L 

+ (flW 9up + 9^(3 9va) D 4 / 2 - 6 (g a/3 <9 M d v + d a dp) D 3 f 2 

- (9au <9 M dp + g a ^ d v dp + gp v d a + g^ d v d a ) D 3 f 2 . (34) 

Note that only the function f 2 , defined in Eq. (J4"|). appears here, in contrast to the 
two-dimensional result. The energy density correlation function in four dimensions 
is given by 

C(x, x') = Cttet? = l — A (V 2 ) 2 D 2 f 2 = l — V 2 V' 2 □ U'f 2 , (35) 

v , ) 7680tt 4 v ; J 7680tt 4 J ' v ; 

where V 2 = □ + <9 2 is the three-dimensional Laplacian operator. This form may be 
used to compute the mean squared average energy density over a spacetime region 
defined by a sampling function F(x). If we define 



p= I d 4 xF(x) :T tt :, (36) 

then 



C= (p 2 ) = J d A xF{x) J d A x' F{x')C{x,x') . (37) 
After an integration by parts, this may be expressed as 

C = -7^4 / d 4 x V 2 DF(x) j d 4 x' V' 2 U'F{x') f 2 (x - x') . (38) 

At first sight, the process of obtaining finite spacetime averages of the correlation 
function may seem mysterious. We start with an expression for C^ap^x.x') which 
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diverges as (x — x')~ 8 as x' — > x, which seems to be a nonintegrable singularity. Yet we 
nonetheless obtain finite integrals of this expression. The reason that this is possible 
is that although C liva p{x : x') is singular as a function, it is a well-defined distribu- 
tion. This is shown by the existence of the expression, Eq. where C fll/a p(x,x') 
is expressed as a sum of derivatives of a function with no more than logarithmic 
singularities. An alternative treatment of the singularities of stress tensor correlation 
functions was given in Ref . 25] . There dimensional regularization was used to render 
the correlation functions finite. In the limit in which n — > 4, where n is the spacetime 
dimension, time-ordered stress tensor correlation functions possess a pole term, which 
can be absorbed in a renormalization involving R 2 and R^ U R^ U counterterms in the 
gravitational action. However, the correlation functions without time ordering, such 
as C^ ua i3(x, x'), have no pole term and are hence finite in dimensional regularization in 
the n — > 4 limit. This is another way to understand why C^ap^x, x') is a well-defined 
distribution, and why the integration by parts method yields finite results. 



B. Averaging over Space and Time - 4D 



Here we will perform a calculation analogous to that in Sect. Ill Bl except involving 
averaging over space and time in four dimensions. The energy density correlation 
function, Eq. may be expressed as 

where r = t — if and r = |x — x'|. As before, we use Lorentzian sampling functions 
of width a in t and in if . The time-averaged correlation function is 

/OO /*00 /*00 

dtg L (a,t) dt'g L {a,t')C(x,x')= I dr g L (a,r) C(x,x') , (40) 
-OO J — OO J — OO 

where a = 2 a. The integrand in the r integral has first order poles at r = ±ia and 
sixth order poles at r = ±r. The integral may be performed by contour integration 
in a way analogous to the integral in Eq. (JTfiJl . The result is 

. _ (3r 2 -a 2 )(r 2 -3a 2 ) 

T 27r 4 (r 2 + a 2 ) 6 ' 1 ' 

Next we wish to average Ct over the spatial directions. Here it will be convenient 
to use a Gaussian sampling function 

gG (l3,x) = -^-e-* 2 /e\ (42) 

in each of the Cartesian space coordinates, x, y, z, x', y', z', and define the spacetime 
average as 

/*oo /*oo /»oo 

C = dxg G (/3,x) / dyg G (p,y) / dzg G ((3,z) 
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x / dx'g G ((3,x') / dy'g G (P,y') / dz' g G {(3, z') C T . (43) 

J— oo J — oo J —oo 

We may use the fact that 

/oo />oo /*oo 

dxg G (P,x) dx'g G (P,x')f(x-x')= / dAx a G (6, Ax) /(Ax) , (44) 
-oo J —oo J —oo 

where Ax = x — x 1 and b = v2/3. This leads to 

/oo roc poo 

dAxg G (b,Ax) / dAyg G (b,Ay) / rfA^ p G (6, Az) C T 
-OO J —OO J —OO 

4 /"°° 

= / drr 2 e- r2 /P 2 C Tl (45) 

where r 2 = (Ax) 2 + (Ay) 2 + (Az) 2 . If we use Eq. (pfTj) . then we can write the spacetime 
averaged correlation function as 

C - r dr r 2 ^ 2 -^)(^-3 a 2 ) m 

The integral in the above expression may evaluated in terms of the error function, 
erf, as 

C = i— (^F fl - erf ( f)l e° 2/V (15b 6 + 90a 2 6 4 + 60a 4 6 2 + 8a 6 ) 

- 2a6(36 2 + 2a 2 )(116 2 + 2a 2 )} . (47) 

Now we wish to discuss the limits in which one sampling length scale is small 
compared to the other. First consider the case of a small spatial scale, 6 C a. The 
exponential factor in Eq. (|46|) guarantees that only values of r < b contribute. Thus 
we can assume that r <C a in the integrand and write 

(3r 2 -a 2 )(r 2 -3a 2 ) ^ 3_ 

(r 2 + a 2 ) 6 ~ a 8 ' 1 ' 

Then we have 

when & < a. This shows that only temporal sampling is necessary in order for C to 
be finite. Equation (f4HJ) may also be derived from the explicit form, Eq. (jUJ), by use 
of the asymptotic form of the error function for large argument. 

Next we consider the opposite limit, where a ^ b. However, C — > oo as a — > 
for fixed, nonzero b. This may be seen from the integral, Eq. (j4fij) . which becomes 
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proportional to J °° drr 6 e r as a — * 0. Alternatively, we can expand Eq. (|4"7|) 
for small a and show that 

C ~ 7/9 1 , 7 , as a — > . (50) 

Thus in four dimensions, averaging over space alone is not sufficient to lead to a finite 
mean squared energy density. This result was obtained previously by Guth [27J and 
byRoura@. 



C. Sampling along a Worldline - 4D 

In the previous subsection, we found that it is possible to take the limit of a 
vanishing spatial sampling scale, so that one is sampling along a worldline. Here 
we will consider that limit for displaced temporal sampling functions. First consider 
Lorentzian sampling functions and let 

/*oo /*oo /*oo 

C(t 0l r) = / dtg L (a,t+t ) / dt' g L (a,t') C(r,r) = / dr g L (a,r-t ) C(r,r) , 



(51) 

where C(r, r) is given by Eq. (jSHJ), and a = 2a. If we were to sample in space with a 
function whose width is small compared to a, the result is the same as setting r = 
in the above expression. More precisely, we perform the integral for nonzero r, using 
the same method as used to obtain Eq. (|lfi|) . and then take the r — > limit. The 
result is 

C(t °' 0) - ^1 + a*y • (52) 

This function has a form similar to that illustrated in Fig. El except that it has three 
maxima and two minima. It is somewhat difficult to graph because the relative sizes 
of the extrema decrease very rapidly. 

In the limit that r = 0, we may write the four-dimensional correlation function as 

C(t, t') = — = l — d? df, \n[(t - i! fie\ . (53) 

v ' ; 2tt 4 (t - t') s 6720vr 4 * * LV ; 1 J K ' 

We can sample the energy density with arbitrary displaced sampling functions and 
write 



DC 



(S(t )S(0)) = dt dt'g(t-t )g(t')C(t,t') (54) 

J — oo J —oo 

-| POO POO 

= "6720^ J „ dt J „ dt> [9t 9{t ~ to)] [d « 9{t ' )] ln [ At2 /^]- (55) 

It should be noted that here we did not use the form of the energy density correlation 
function, Eq. (|35j). which follows from the covariant form. Instead, we let r — > 0, 
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FIG. 5: The graph of K(to) in four dimensions as a function of for the compactly 
supported sampling function given by Eq. I|56j) . in units where a = 1. 



and then expressed the result in terms of time derivatives of a logarithmic function. 
A more rigorous approach would be to average Eq. over both space and time, 
and then let the widths of the spatial sampling functions go to zero. However, this is 
difficult to do explicitly with general sampling functions. The equivalence of the two 
approaches needs to be studied more carefully. 

Let us next consider a compactly supported sampling function given by 



9(t) = ^ (t - a/2) 4 (t + a/2)* 



(56) 



for \t\ < a/2, and g(t) = for \t\ > a/2. Note that g(t) and its first three derivatives 
vanish at t = ±|a, so all surface terms vanish when we integrated by parts in Eq. (fo3j) 
to obtain the second form for (S(t )S(0)). We may again define K(t ) by Eq. (J2HJ) 
and evaluate it numerically. The result is plotted in Fig. 

As in two dimensions, there are regions of correlation and of anticorrelation as to 
increases. However, the behavior in four dimensions is more complicated, with three 
maxima and two minima. This appears to be due to the greater number of derivatives 
of the sampling function in Eq. (|53j) . as compared to Eq. (J23j) . 

In Appendix O we show that in two and four dimensions 



and that 



poo 

/ K(to)dt = , 

<s 2 (o)>>o, 



(57) 
(58) 



14 



for a general g(t). From Eq. we can also explicitly verify that / °° K(t Q )dt = 
for the Lorentzian sampling function. 



IV. SUMMARY 

In this paper, we have presented covariant expressions for the Minkowski vacuum 
stress tensor correlation function in two dimensions, Eq. (J2J), and in four dimensions, 
Eq. (|34|). These expressions are of the form of a sum of terms, each of which is a 
derivative of a scalar function with logarithmic singularities in the coincidence limit. 
These expressions allow one to write spacetime averages of the correlation function 
as finite integrals. We explicitly evaluated such averages of the energy density in 
two dimensions using Lorentzian sampling functions in both space and time. The 
resulting expression, Eq. ()17j) . is symmetric in the spatial and temporal sampling 
widths, and is finite as either width goes to zero with the other width fixed at a 
nonzero value. 

We next studied the correlations of the sampled 2D energy density along a world- 
line using displaced sampling functions. This reveals the correlation and anticorre- 
lation of measurements of the energy density in overlapping intervals. The result is 
illustrated in Fig. H]for a compactly supported sampling function. When the intervals 
nearly overlap, the two measurements are positively correlated, as expected. When 
the overlap has decreased somewhat, the two measurements become anticorrelated. 
This can be interpreted as telling us that if we find energy density of one sign on 
the first measurement, we should find the opposite sign on the next measurement. 
Finally, as the intervals become disjoint, the measurements are again positively cor- 
related. Furthermore, we show that for an arbitrary sampling function, the net area 
under the correlation graph, e.g., the one depicted in Fig. HI is equal to zero. It is 
hoped that further investigation will elucidate this interesting behavior. 

The analogous calculation in four dimensions yields similar results. However, in 
this case there are two regions of anticorrelation and three of positive correlation. 
The fluctuations in the averaged energy density remain finite in the limit that the 
spatial width vanishes, but not in the limit that the temporal width goes to zero. 
Thus in four dimensions, the averaged energy density correlation function requires 
averaging in time to be finite. 

There is a vaguely analogous result concerning quantum inequalities on the av- 
eraged expectation value of the stress tensor in an arbitrary state. There are finite 
lower bounds on the expectation value of the energy density averaged on a worldline 
in both 2D and 4D, and on the spatial average in 2D. However, the spatial aver- 
age in 4D has no lower bound [29]. The search for a deeper link between quantum 
inequalities and the vacuum stress tensor correlation function is a topic for future 
research. 

Another question which needs to be explored further is that of the physical effects 
of the passive metric fluctuations driven by vacuum stress tensor fluctuations. One 
approach is that adopted in Ref. [2O] where the Raychaudhuri equation was used as 
a Langevin equation to study the luminosity fluctuations and angular blurring of a 
distant source produced by passive metric fluctuations. The case of the Minkowski 



15 



vacuum was briefly discussed in Ref. [2(|, where it was found that the natural quan- 
tum uncertainty in the test particles used to probe the fluctuating geometry tends to 
hide the effects of the metric fluctuations. However, this does not necessarily mean 
that these fluctuations are in principle unobservable. This is another question for 
further study. 
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APPENDIX A 

In this appendix, we give the derivation of the stress tensor correlation function in 
two-dimensional spacetime in the Minkowski vacuum state. We first start with the 
form of the stress tensor for a massless, minimally coupled scalar field: 

T,a, = ^,n<j), v ~ -xg^u ,P 0, P • (Al) 
From this expression we find the correlation function 

C^ a p(x,x') = (: T^(x) : : T a/3 (x') :) = (: d^(j)d v (j) : : d a ,(j)dp,(j) :) 

~ 2^ ( : d p (pdp4> : : d a ,(f)dp><f) :) - -g a p (: d^d^ : : 0(9^0 :) 
- \g^9ap (■ d p <j)d p( f) : : d°' d a , :) . (A2) 

Here unprimed indices refer to the point x and primed indices to x' . Next we use the 
identity 

(: 01 2 ■ ■ 03 04 ■) = (01 03) (02 04) + (01 4>A) (02 03) , (A3) 

where the 0j are quantum fields or derivatives of quantum fields. From this identity, 
we can show that 

(: ^0 d v (j> : : ^0^0 :) = (3„ d a , D){d v dp, D) + (3„ dp, D){d v d a , D) , (A4) 
where 

D = D(x, x') = (0(x)0(x')) (A5) 

is the two-point function. We can express the correlation function in terms of deriva- 
tives of D as 

Cfvap (x, x) = (<9 M d a > D) (d u dp, D) + (d^ dp, D) (d u d a , D) 

~ 9lxv {d" d a , D) [d p dp,D)- g a p (d, d°' D) [d v dfj, D) 

- \g^g a p(ded°' D)(d p d a ,D). (A6) 
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An equivalent expression for the case of a massive, nonminimal scalar field has been 
given by Martin and Verdaguer. (See Eq. 3.42 in Ref. [H3].) The analogous expression 
for the electromagnetic field is given in Ref. jl8l ]. 

Up to this point, our treatment applies to spacetimes of any dimensionality. Now 
we specialize to two-dimensional Minkowski spacetime. There is an infrared diver- 
gence in the two-point function for a massless scalar field in the Minkowski vacuum 
state in two dimensions. Thus, the field must either have a nonzero mass, or else the 
only physically allowed states are ones which break Lorentz invariance [30]. Fortu- 
nately, the details of either approach have no effect on our results. If we let the scalar 
field have a small mass m, then the two-point function is given by 

D = — ln(cm 2 Ax 2 ) (A7) 

47T 

in the limit that m~ 2 ^> Ax 2 . Here c is a dimensionless constant and Ax 2 = (x M — 
x m )(x /Ji — x'p). Because the stress tensor correlation function depends only upon 
derivatives of D, it is independent of c and m. The second derivative of D is 

where Ax M = x^ — x' . We can now combine Eqs. (jA6|) and (jA8|) to obtain an explicit 
expression for the stress tensor correlation function in two dimensions: 

1 r 8 

C^apix,^) = — ^ ~ Ax M Ax„Ax a AX(3 



4vr 2 
2 



(Ax 2 

(g m Ax v Axp + Ax v Ax a + g va Ax^ Axp + g vP Ax M Ax c 



(Ax 2 ) 3 

+ ( Ax2 y (sw* 9up + g^ gua - g^u g a p) (A9) 

We next wish to express C fJiUa p(x,x') as a sum of derivatives of scalar functions. 
Lorentz symmetry suggests that these be functions of Ax 2 . Let / = /(Ax 2 ). Then 
the derivatives of / are 

d ll f = 2Ax ll f, (A10) 
d, d v f = 2 g, v f + 4Ax /t Ax, /" , (All) 
<9 M d v d a f = 4(g^ Ax a + g m Ax u + g va Ax^f" + 8Ax^ Ax v Ax a f" , (A12) 

and 

d u d a dpf = A(g^ a g vfi + g^ g ua + g^ g aj3 )f" + 8(g^ a Ax v Ax p + g^ Ax v Ax a 
+ g va Ax M Axp + g uf3 Ax M Ax a + g^ Ax a Axp + g aj3 Ax^ Ax u ) f" 
+ 16A^ Ax v Ax a Ax p f"" . (A13) 

Here primes denote derivatives of / with respect to its argument. We will also need 
some expressions involving the wave operator: 

□/ = d^f = 2nf + 4Ax 2 f , (A14) 
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□□/ = 4n(n + 2)f" + lQ(n + 2)Ax 2 f" + 16(Ax 2 ) 2 /"" , (A15) 

and 

d, d u □/ = 4(n + 2)g, u f" + 8[^ Ax 2 + (n + 4) Ax M Ax,] /"' 

+ 16Ax M Ax v Ax 2 f" , (A16) 

where n is the dimension of the spacetime. 

Our goal is to express C /JiVa p(x,x / ) as a sum of derivatives acting on one or more 
choices of /. Because d^f\ = — dn? /i, we can write our results in several equivalent 
forms, but here and in Appendix[B]we will use derivatives with unprimed indices. If / 
is dimensionless, then in two dimensions we will need four derivatives in each term in 
order that C flllo ,p(x, x') has dimensions of length -4 . There are five fourth-rank tensors 
that we can form which have the correct dimensions and symmetry properties: 

d li d u d a d /3 f, (A17) 

(gapdpd,, + gp,d a dp)Uf , (A18) 

(9au dp dp + Qafi d v dp + g Pv <9 M d a + d v d a )□/ , (A19) 

g^g^DDf, (A20) 

and 

(9^a9u/3 +sy 1( sflwpD/. (A21) 

We would like / to have an integrable singularity at Ax 2 = 0, so a natural choice 
is a power of a logarithmic function. First consider 

h = HAx 2 /i 2 ) , (A22) 

where I is an arbitrary constant with dimensions of length. However, Ofx = in two 
dimensions, so the only nonzero tensor from the above list which can be formed from 
fi is 

96 

d^d u d a d p fx = ~ ( Aa . 2 )4 Ax u Ax a Axf3 
16 

+ ^ x2 y (9iu* Ax » Ax P + 9nP Ax » Ax * + 9ua Ax M Ax^ 

+ g v p Ax M Ax a + g^ Ax a Ax p + g a/3 Ax M Ax v ) 
4 

_ (A~2Y2 (9^a 9u/3 + 9^/3 9ua + 9/j.u 9afi) ■ (A23) 

This is not sufficient to form C Mi , aj g(x, x'), so we need another choice of /, which we 
take to be 

f 2 = ln 2 (Ax 2 /£ 2 ) . (A24) 
From Eqs. (|A14|) and (|A15j) with n = 2, we find 

□/ 2 = ^ (A25) 
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and 

32 

nn ^=(A^r < A26 > 

This allows us to form four tensors from f 2 with the correct symmetry properties and 
dimension: 

g af3 DDf 2 = 32 g af3 , (A27) 

(sVa fl^ + flr^s 5W) EC/a = 32 (3^ ^ + g^ g ua ) ^ 2 y , (A28) 

Cfc*3 M ft, +^9 a a^p/ a = -32 ^§ + 64 ^ AX ° Ax ^f Ag " AX - , (A29) 
and 

g m gup + gup g va 



{g<xv dp dp + g ail d v dp + gp v dp d a + gpp d v d a )Uf 2 = -32 



(Ax 



2\2 



64 

+ ^ Aj;2 ys \9na + Ax v Ax a + g ua Axp Axp + g vP Axp Ax a ) . (A30) 

Note that <9 M d v d a dp f 2 is not a suitable term because it contains logarithmic pieces 
that do not appear in Cp Ua p(x, x') and which cannot be cancelled by any other terms. 
This leaves us with five tensors from which to form the stress tensor correlation 
function. 
Let 

C^ap {x,x') = 3g ^ 2 ci 8^ d v d a dp fx + c 2 fiv nOf 2 

+ c 3 g^g + g^p g ua ) UUf 2 + c 4 {g a p dp + g^ d a dp)Uf 2 

+ c 5 (g av dp dp + g a p d v dp + gp v dp d a + gpp d u d a ) Df 2 . (A31) 

If we insert Eq. (IX23|l and Eqs. (IX27|) - (IX30|) into this expression and compare with 
Eq. (|A9J) . we find five conditions on the five coefficients. The unique solution of these 
conditions gives 

c\ — — 8 , c 2 = — C4 = —2 , and C3 = — c 5 = 1 . (A32) 

As a check, the correlation function may be shown explicitly to satisfy the conserva- 
tion law 

& Cp va p{x, x') = d a ' Cp Ua p{x, x') = . (A33) 



APPENDIX B 

Here we repeat the derivation in the previous appendix for the case of four- 
dimensional Minkowski spacetime. The general form, Eq. (|A6|) . for the correlation 
function still holds, but the two-point function for a massless scalar field is now 

D = ra- (B1) 



19 



If we insert this form into Eq. (jA6|) . we find the four- dimensional analog of 
Eq. 



32 

C^ a p(x,x f ) = -^—^ Ax^ Ax u Ax a Axp 



4tt 4 
4 



(g^ a Ax u Axp + g^p Ax v Ax a + g va Ax M Axp + g v p Ax^ Ax a ) 



(Ax 2 ) 5 
g 

{9 iiv &x a Ax p + g a p Ax^ Ax 



+ ^ x 2\a (5W 9up + g^ gua + 4sv g a p) (B2) 

In four dimensions, the correlation function has dimensions of 1 /length 8 . Thus 
any expression involving derivatives on a dimensionless function will require eight 
derivatives. Because there are only four free indices, there will have to be at least 
two wave operators. This eliminates the logarithm function fx, Eq. (jA22|) . because 
in four dimensions 

□□/i = 0. (B3) 

However, the squared logarithm function f 2 may be used to form the following five 
tensors with the correct dimensions and symmetry: 

d„d v d a d p nnf 2 , (B4) 

{g a p dfj. d u + g^ d a dp)D 3 f 2 , (B5) 

(flW <9 M dp + g afl d v dp + gp v <9 M d a + g PfJL d v d a )D 3 f 2 , (B6) 

g^u g a p D 4 f 2 , (B7) 

and 

(9ica 9vp + QnP 9va)0 4 f 2 . (B8) 

We may repeatedly use Eqs. (jA14|) and (jA15|) with n = 4 to demonstrate that, in 
four-dimensions, 

32 

DD/2 = -(A^' < B9) 

D3 *=-(Jfr < B10 > 

and 

From these expressions, we may show that 

61440 

df, d u d a dp UUf 2 = ~ ^ x 2\e Ax f Ax u Ax a Axp 
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6144 

+ ^ g2 )5 ( 9 ^ a AXu AX P + Brf AXu AXa + 9ua AX P 

+ g v p Ax M Ax a + g^ v Ax a Axp + g a p Ax M Ax u ) 
768 

~~ (Ax 2 ) 4 a 9 1/13 ^ ^ ® ua ^ ^ ® ' 



(B12) 



(fiW df, dp + g afJi d v dp + g pv d^ d a + d v d a )D 3 f 2 



3072 



( V* g v p + g»p g 



12288 



(Ax 



2\5 



(g m Ax v Ax/3 + g^p Ax u Ax 

a ~r gua 

Ax^ Axp + g v p Ax M A 



(B13) 



and 



(gap <9 M d v + g^ d a d p )n*h = 3072 



gjw gap _ g^u Ax a Axp + g a p Ax M Ax v 
(Ax" ' 



2^4 



(Ax 2 ) 5 



We now express the correlation function as a sum of the tensors formed from f 2 



as 



1 



61440 vr 4 



ci df, d v d a dp UUf 2 + c 2 g^ g a p 4 f 2 



+ c 3 {g m 9vp + 5V/3 Qua) D 4 /2 + c 4 {g a p <9 M d u + ^ d a dp)D 3 f 2 

+ c 5 (# ai , d^ dp + g afl d v dp + gp v d^ d a + g Pfl d v d a ) D 3 f 2 . (B15) 



If we insert the explicit forms for these tensors and compare with Eq. (jB2|) . we again 
find five conditions on the five coefficients, leading to the solution 

C\ = —8, c 2 = — c 4 = —6, and c 3 = — c 5 = —1 . (B16) 

As required, the correlation function has a vanishing divergence on any index. 



APPENDIX C 

In this appendix, we will prove Eqs. ()32|) and (}3*3*j) in both two and four dimensions. 
We will proceed by first showing that 



K(t )dt = 0. (CI) 

> 

Then we will prove that If (to) is a symmetric function, and hence show that 

K(t )dt = 0, (C2) 
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as well. 

Let g(t) be an arbitrary smooth sampling function. From Eq. (J2%j) or (J51j) and 
Eq. pSjl. we have 

I /"OO /*oo 

K{t ° ) = WmJ dt9{t ~ to) J dt '9(t')C(t-t'). (C3) 



Then 



/oo i /-OO /-oo /*oo 

K(t )dt = — — / / dt g(t-t ) / dt' g{t') C{t-t') 

-I /»oo /»oo 

= wmL dt ' 9{t ' ) J-J tc{t - t ' h (C4) 

where we have interchanged the order of integrations, and used the fact that for 

y = t-ta, 



dt g(t-t ) = - dyg(y)= dy g(y) = 1 . (C5) 

) ./ oo •/ — oo 

However, if we can write C(t — t') = dF(t — t')/dt, where F(t — t') — > as t — > ±cxd, 
then 



which in turn implies that 



/oo 
dt <7(f - = [F(t - t')t + Z = , (C6) 
-oo 

/oo 
K(t )dt = 0. (C7) 
-oo 

Recall that in two dimensions, the worldline vacuum correlation function is C(t— t') = 
l/[47T 2 (t - t') 4 ], and in four dimensions it is C(t - t') = 3/[2vr 4 (t - t') 8 ], so in both 



cases the condition Eq. (jC6|) is satisfied. Note that in four dimensions, it is necessary 
to assume that we set the spatial separation r in Eq. (|39j) to zero and then average 
over time, as discussed in Sect. IIII CI 

We now show that K(t ) = K(—t ). Let us write 

/oo /»oo 
dtg(t-t ) / dt'g(t')C(t~t') 
-oo J — oo 

OO POO 

dt g{t) I dt'g(t') C(i + t - t') 



oo J — oo 

OO POO 



/oo poo 
dtg{t) I dFg(V + t )C(t-F), (Ci 
-oo J — oo 



where we have let t = t — 1 , so t = t + to, and t' = t' — 1 . If we now let t' — > t, t — > t', 
we have 



(5- 2 (0))A:(to) = / dtg(t + t ) / dt' g{t')C{t' -t) 

J —oo J —oo 

= (S 2 (0))K(-t ), (C9) 
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where we have used the fact Cii! — t) = C(t — £'). Note that the symmetry of K(to) 
depends only on that of C and does not assume that the sampling function g(t) is 
symmetric. Thus since 



K(t Q )dt Q = 0, (CIO) 
and K(t ) is symmetric, it also follows that 

POO 

/ K(t ) dt = . (Cll) 
Jo 

In order to determine whether a fluctuation is correlated or anti-correlated with 
itself, we must determine the sign of (S 2 (0)) in the general case. We would expect that 
a fluctuation should be correlated with itself, and thus that (5 2 (0)) > 0. This can be 
proven from the fact that 5(0), as defined by Eq. (J22|) . is a self-adjoint operator pul ]. 
Let \ip) be the state under consideration, which in our case is the Minkowski vacuum. 
Then 

|*)=£(0M (C12) 
is a well defined state vector with positive norm. Thus we have 

H^H 2 = (^15^0) 5(0) |V) = (5 2 (0)) > 0. (C13) 
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